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1. Introduction 

Quantum chromodynamics (QCD) describes the interaction between quarks and gluons which 
is responsible for the strong nuclear force. Recent advances in QCD using computational lattice 
methods have made it possible to accurately predict the spectrum and properties of many isolated 
hadrons. Unfortunately, lattice QCD calculations of nuclear and neutron matter or even few-body 
systems beyond two nucleons are presently not possible. The most significant challenge in such 
simulations would be to overcome the exponentially small signal-to-noise ratio caused by the sign 
and complex phase oscillations for simulations at large quark number. 

Nuclear lattice simulations based on EFT provide an alternative method to describe few- and 
many-body systems at low energy without losing connection to QCD. The lattice EFT approach 
addresses the few- and many-body problem in nuclear physics by applying non-perturbative lattice 
methods to low-energy nucleons and pions. The effective Lagrangian is formulated on a spacetime 
lattice and the path integral is evaluated by Monte Carlo sampling. Pions and nucleons are treated 
as point-like particles on the lattice sites. By using hadronic degrees of freedom and concentrating 
on low-energy physics, it is possible to probe large volumes and greater number of nucleons than in 
lattice QCD. After a brief overview what has been done in this field so far we present some results 
of our recent studies of the two nucleon system [Q] and neutron matter at subleading order. 
Accurate description of two-nucleon phase-shifts and ground state energy ratio of dilute neutron 
matter up to corrections of higher orders show that lattice EFT is a promising tool for quantitative 
studies of low-energy few- and many -body systems. 



2. Lattice EFT: previous achievements 

Lattice EFT is a rather new and fast developing field. Here we give a brief overview on what 
has been done in this field so far. For a comprehensive discussion the reader is referred to [Q]. The 
first lattice study of nuclear matter was carried out in the early nineties by Brockman and Frank [^] 
using a momentum lattice and based on the hadrodynamics model of Walecka [^. The first nuclear 
lattice calculation based on EFT was carried out by Miiller et al. [^. They looked at infinite nuclear 
and neutron matter at nonzero density and temperature. Later a series of analytical studies were 
carried out: Chen and Kaplan ^ showed the absence of sign oscillation for nonzero chemical 
potential in the Hubbard model. Non-linear realization of chiral symmetry with static nucleons on 
the lattice was discussed by Chandrasekharan et al. [^. Also ChPT within the lattice regularization 



was considered by several groups [|10|, [11|, |12|]. This was followed by the first many-body lattice 
calculation using chiral EFT [13]. Since that time a number of lattice calculations for cold atoms 
and low-energy nuclear physics were carried out. See [Q] for a review article. It is important to 
note that in the low-energy sector the phase region accessible by lattice EFT is much broader than 
in lattice QCD. Severe sign oscillation problem limits the accessibility of finite density lattice QCD 
simulations. In contrast, sign oscillations in nuclear lattice EFT are strongly suppressed due to the 
approximate SU(4)-symmetry in the two-nucleon sector. One can show explicitely that SU(4)- 
symmetric nuclear EFT does not have a sign problem. SU(4) symmetry breaking leads to small 
sign oscillations which, however, turns out to be not severe. 
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NLO 


N^LO 


N^LO 


Exp 


Ed [MeV] 
As [fm-V2] 


-2. 171. ..-2. 186 
0.0256... 0.0257 
0.868... 0.873 


-2.189... -2.202 
0.0255... 0.0256 
0.874... 0.879 


-2.216... -2.223 
0.0254... 0.0255 
0.882. ..0.883 


-2.224575(9) 
0.0256(4) 
0.8846(9) 



Table 1: Deuteron properties at NLO, N^LO and N^LO compared to the data. Here, is the binding 
energy, Tjd the asymptotic D/S ratio and the strength of the asymptotic S-wave normalization. The data 
for are from [Q, for rj^ from 1^] and for A5 from [p5[]. 



3. Nuclear EFT 

Let us now give a brief introduction to the basic foundations of our approach. The low-energy 
properties of hadronic systems are, in principle, accessible in lattice QCD. This method is, however, 
very expensive, especially if one wants to consider few- and many-nucleon systems. Alternatively, 
we can exploit the spontaneously broken approximate chiral symmetry of QCD which implies 
the existence of light weakly interacting Goldstone bosons. In the SU(2) sector, we identify the 
Goldstone bosons with pions. Since the interaction between the Goldstone bosons is weak one 
can apply perturbation theory, where the expansion parameter is not a coupling constant but small 
momenta and masses of the Goldstone bosons divided by the chiral symmetry breaking scale A^,. 



This systematic procedure is called chiral perturbation theory [14] and reproduces (as explicitely 



proved by Leutwyler 0150 in mesonic sector) order by order the original QCD Green-functions. 

ChPT has been extended to one nucleon sector. In the two and more nucleon sector additional 
problems appear. Due to the existence of nuclear bound states, the strict perturbative procedure 
breaks down. As shown by Weinberg [|l6|], the power counting is violated by nucleon-nucleon 
(NN) cuts. He suggested to construct perturbatively a so-called chiral effective potential which, 
per construction, excludes the NN cuts and, for this reason does not violate the power counting 
of ChPT. To describe nuclear observables in the two-, three- or more-nucleon sectors one should 
numerically solve the Lippmann-Schwinger, Faddeev or Faddeev-Jakubowsky equations, respec- 
tively, with the chiral effective potential as an input. 

Chiral effective potential has been extensively studied in the last decade up to next-to-next-to- 
next-to-leading order (N^LO) in chiral expansion (for extensive discussion see [17]). At this order 



two leading-order (LO), seven subleading order (NLO) and fifteen N-^LO unknown low energy 
constants have been fitted to low energy nucleon data [|l^, [l^, At this order in the chiral 
expansion, one observes an accurate description of all NN low-energy observables, see Figs. [I], ^ 
and Table |. 

4. Nuclear EFT on the lattice 

Once the chiral nuclear forces are determined and the low energy constants appearing in the 
nuclear forces are fitted (in the two and thiee-nucleon sector) one can make predictions in the 
four- and more-nucleon sectors based on chiral EFT. However, explicit numerical treatment of the 
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50 100 150 200 250 50 100 150 200 250 50 100 150 200 250 



Lab. Energy [MeV] Lab. Energy [MeV] Lab. Energy [MeV] 

Figure 1: S-, P- and D-waves np phase shifts. The dashed, Hght shaded and dark shaded bands show 
the NLO, N^LO and N-'LO [jl|] resuhs, respectively. The dashed Hne is the N^LO result of Ref. [Q. The 
filled circles (open triangles) depict the results from the Nijmegen multi-energy PWA pc[ , pi] ] (Virginia 
Tech single-energy PWA p^). 



Jakubowsky equations for more than four nucleons is a very difficult task. To solve the many- 
body problem we propose to put the chiral effective potential on the lattice and apply the powerful 
Monte-Carlo techniques which are already developed to high degree. In this framework, nucleons 
are represented as point-like Grassman-fields and pions as point-like instantaneous (in order to 
reproduce the chiral potential) pseudoscalar fields. Typically, our calculations are carried out using 
the lattice length L ~ 20 fm and the lattice spacing a ~ 2 fm which corresponds to the cutoff 
A = tt/a: ~ 300 MeV. The correlation function for A nucleons in the Euclidean space is defined by 

ZA{t) = {^AWxp{-tH)\^'A), (4.1) 

where the states \^a) refer to the slater determinants for A free nucleons, H is the Hamiltonian 
of the system and t the Euclidean time. The ground state energy of the A-nucleon system can be 
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Figure 2: np differential cross section and vector analyzing power at E\ah = 50 MeV. The Nijmegen PWA 
result is taken from [pT[]. 



derived from the asymptotic behavior of the correlation function for large t. 

£0 = -lim— lnZA(0- (4.2) 
dt 

Expectation value of a normal ordered operator t& can be derived in a similar way by 

{^'a\^\^a) = lim Zf{t) = {^>A\exp{-tH/2)i^exp{-tH/2)\^'A >, (4.3) 

where the states denote the ground states of A-nucleons system. It is convenient to de- 
scribe NN contact interactions by standard bilinear nucleon density operators using the Hubbard- 
Stratonovich transformation. Using the relation 

exp{p^/2) J dsexp{-s'^/2-sp) (4.4) 

one can express terms quadratic in the nucleon density operator p as terms linear in p in the 
presence of auxiliary background fields. In this representation, the full correlation function is 
related to the path integral over pions and auxiliary fields, 

){WA\M^'"-'\7i,,st)---M^^\7ij,st)\^'A). (4.5) 

7=1,2,3 i 

Here Sjt^i and 5„ are free actions for pions and auxiliary fields Si and M^"^ is a transfer matrix 
defined as an n'th step in the temporal direction. Note since we only have linear nucleon density 
operators in the action the amplitude 

{^>a\m'^'^'-'Hki,si)---M^^\ki,si)\^>a) (4.6) 

is just a slater determinant of single nucleon matrix elements ^ y with i,j= 1, . . . ,A. 



5. Lattice EFT at leading order 

To be specific, we give here the leading order action starting with the free theory. The presen- 
tation here is somewhat sketchy. For an extensive discussion see [26]. The free auxiliary fields and 
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pion actions are given by 

Sss{s,si) = ^i:.(n)' + ^ II^/(n)^ S,,{7ti) = ^ £ ^;rKn)(-A + M2);r,(n), (5.1) 

n 



'1=1 n 



1=1 n 



where M„ is the physical pion mass, / denotes isospin indices and (Xt = at/a. For nucleons we use 
O(a^) improved free lattice Hamiltonian defined by 



1 



"■Sl'^Sl'-lJ 



(5.2) 



where m is the nucleon mass, the operators aj (n^) and aij{ns) are the nucleon creation and anni- 
hilation operators, Hg are spatial coordinates, are spatial unit vectors, the indices i and j stay for 
spin and isospin indices, respectively, and the coefficients fk read: 



_ 49 _3 ^ 1 

/o,i,2,3-y>-4>^>-Y^- 



(5.3) 



To define the interactions we introduce the nucleon-density operators with different spin/isospin 
polarizations 

p°^«(^,) = £4/n,)a;,y(n,), pf'"{n,) = £ «;,-,(n,)[T/]/ja,-j("^)> (5-4) 



Pif^ns)= £ 4,/("^)[<^5],-',/[T/]/ja;j(n,). (5.5) 
The transfer matrix for n^-th step has, besides the free part, two important contributions: 

= :exp j-77freea^-^Ei:Vs;r,(n„nOP|,?(«.) 

L ^ S,I Hs 



(5.6) 



Here :: denotes normal ordering. The first long-range contribution includes the instantaneous pion- 

nucleon interaction and describes the one-pion-exchange in the leading-order effective potential. 
The second short-range contribution corresponds to the NN contact interactions. The low-energy 
constants C and C/ fitted to Nijmegen PWA appear to have different signs: 



C < 0, C/ > 0. 



(5.7) 



With these signs the pion-less theory appears to have no sign-oscillations if the number of protons 
and neutrons are equal and they stay pair-wise in isospin-singlet states. In this case the multipli- 
cation with T2 of the single-nucleon matrix elements ^ from left and right is well defined and 

gives 

T2^T2 = .^*. (5.8) 

For this reason, the determinant of ^ appears to be real: 



(5.9) 
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Since T2 is antisymmetric, the eigenvalues of ^ are doubly degenerate. This leads to a positive 
slater determinant [|8, 27] 

det^>0. (5.10) 

The introduction of pions causes small sign-oscillations which, however, are not severe and appear 
to be suppressed. 

To perform our simulations in a most efficient way we, in addition, exploit the approximate 



SU(4)-Wigner [28| symmetry in the NN system. The symmetry transformation is given by inde- 
pendent rotation of spin and isospin degrees of freedom. 

5A^ = a^va'^T^'A^ with o^ = {\,d) and T'' = (l,f). (5.11) 

One can show that in the limit where the NN S-wave scattering lengths approach infinity the two- 
nucleon system becomes invariant under the SU(4)-transformation ||29|]. The SU (4) -breaking cor- 
rections come from the finite scattering length and higher order terms in the chiral expansion: 

SU(4) -breakingterms ~ -3!— - (5.12) 

Since the NN scattering lengths 

a(i5o) = (-23.758 ±0.010)fm = (5.424 ±0.004)fm (5.13) 

are very large, the SU (4) -breaking corrections appear to be small. This fact can be used to improve 
the performance of our lattice simulations. The SU(4) symmetric transfer matrix is given by 



m("') =: exp 



-//tree Or + ^^C^^ ^(H,, «,)p''^ (^ 



(5.14) 



In this case there are no sign-oscillations for even number of nucleons [ ]30| ] and we do have only 
one auxiliary field such that the simulations are much cheaper. Although there is no positivity 
theorem for odd numbers of nucleons, sign oscillations appear to be suppressed also in systems 
with odd number of nucleons because it is only one particle away from an even system with no 
sign-oscillation. Since the final result is closed to the one produced by SU(4)-symmetric simulation 
we divide our simulations in three parts. To simulate an expectation value of some observable we 
use SU(4)-symmetric transfer matrices in the first and the last steps in order to filter the low- 
energy signal and after filtering start the simulation with realistic transfer matrices. A schematic 
overview of the transfer matrix calculation is shown in Fig. 0. 



For our nuclear lattice simulations we use the hybrid Monte-Carlo (HMC) method [|31|]. We 
introduce the conjugate fields pji,, p^, Ps, and use molecular dynamics trajectories to generate new 
configurations for the fields pjc,, ps, Psj, Tij, s, sj which keep the HMC Hamiltonian 

//hmc = ^5 [pl{n) + pUn)] +pUn)^ +V{nj,s,sj), (5.15) 

constant, where the HMC potential is defined by 

V{7i:i,s,si)=Syt^ + S,,-log\det.^\. (5.16) 

Upon completion of each molecular dynamics trajectory, we apply Metropolis accept or reject step 
for the new configuration according to the probability distribution exp (—//hmc)- This process of 
molecular dynamics trajectory and Metropolis step is repeated many times. 
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2Lf + Lf. L, + L,. 
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expectation value 

/ 

full LO O full LO 
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Figure 3: Overview of the various pieces of the transfer matrix calculation. 







r3H[fm] 


£4He[MeV] 


''4He[fm] 


Simulation 
Experiment 


-8.9(2) 
-8.482 


2.27(7) 
1.755(9) 


-21.5(9) 
-28.296 


1.50(14) 
1.673(1) 



Table 2: Experimental and nuclear lattice simulation results for binding energies and root-mean-square 
radius of triton and helium-4. 



6. Leading-order results 

With the presented method we performed nuclear lattice simulations on JUBL/JUGENE su- 
percomputer at Forschungszentrum Jiilich. Already at leading order we get promising results for 



binding energies, radii and density correlations for the deuteron, triton and helium-4 [|26|]. Numer- 
ical results on a 5^ lattice for triton and helium-4 are shown in Table ^ The triton binding energy 
agrees with experiment within 5% and the triton root-mean-square radius is accurate to 30%. The 
binding energy for helium-4 is within 25% of the experimental value while the root-mean-square 
radius agrees within 10%. Our results for the triton nucleon density correlations are shown in 
Fig. ^ We also studied the feasibility of simulations for light nuclei with up to eight nucleons and 
observed that for A < 8 the CPU time scales approximately linear with A. 



7. Next-to-leading-order results 

At NLO there appear 9 low energy constants (LECs) which we fitted to the Nijmegen NN 



scattering data. Elastic scattering phase-shifts on the lattice are related by Liischer's ||32|, |33|, |34| ] 
formula to the energy levels of two-body states in a finite large volume cubic box with periodic 
boundary conditions. While this method is very useful at low momenta, it is not so useful for 
determining phase shifts on the lattice at higher energies and higher orbital angular momenta. Fur- 
thermore, spin-orbit coupling and partial-wave mixing are difficult to measure accurately using 
Liischer's method due to multiple-scattering artifacts produced by the periodic cubic boundary 
conditions. In ||l|] we proposed a more robust approach to measure phase shifts for two nonrela- 
tivistic point particles on the lattice using a spherical wall boundary. The basic idea is to impose 
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G^Jx,y,0) 




X (lattice units) 



Figure 4: The nucleon density correlation for the triton in xy-plane. 



a hard spherical wall boundary on the relative separation between the two interacting particles at 
some chosen radius. The reason for this spherical wall is to remove copies of the two-particle in- 
teractions due to the periodic boundaries on the lattice. This additional boundary condition allows 
for a direct extraction of the phase-shifts and mixing angles from the finite-volume spectrum. For 
more details see [|l|]. 

Using the spherical wall method we determined the values of 9 LECs by matching three S- 
wave, four P-wave scattering data points, as well as deuteron binding energy and quadrupole mo- 
ment. In Fig. H we show NN S-wave phase-shifts and mixing angles for two different actions, called 
LOi and LO2. The action LOi is the one presented in Eq. In the action LO2 the contact inter- 
actions are smeared by a Gaussian. The two actions are identical at leading order and differ only 
by higher-order terms. The main motivation to introduction the Gaussian smearing was to cure a 
multi -particle clustering instability at coarse lattice spacing present in simulations with LOi and to 
estimate a systematic error coming from higher-order corrections see [26] for extended discussion. 
As can be seen from Fig. ||, our lattice simulation results are in a good agreement with the partial 
wave results for momenta smaller than 80 MeV. Deviations between the two results for different 
actions appear merely at larger momenta and are consistent with the expected higher order effects. 



8. Dilute neutron matter 

As a first application at NLO we simulate dilute neutron matter in a periodic box [^]. We 
probe the density range from 2% to 8% of normal nuclear matter density. Neutron-rich matter at 
this density is likely to be present in the inner crust of neutron stars. The Pauli suppression of 
three-body forces in dilute neutron matter makes it a good testing ground for chiral EFT applied to 



9 



Nuclear effective field theory on the lattice 



Hermann Krebs 




-10 
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PcM (MeV) 



-10 
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PCM(MeV) 



Figure 5: NN S-wave phase shifts and mixing angles versus center-of-mass momentum with actions LOi 
and LO2. 



many-nucleon systems. Neutron matter at kf ~ 80 MeV, where 

kF = ^{37i^Ny^^ (8.1) 

is Fermi momentum, is close to the so-called idealized unitary limit. In this limit the S-wave 
scattering length is infinite and the range of the interaction is zero such that the scattering amplitude 
is as strong as possible. At lower densities corrections due to the finite scattering length become 
more important while at higher densities corrections due to effective range start to dominate. In the 
unitary limit the ground state has no dimensionful parameters other than the particle density and so 
the ground state energy of the system should obey a simple relation 

Eo = ^4"' (8-2) 

where is a dimensionless measurable constant. Due to its universal nature, the unitary limit can be 
studied in atomic systems. Ultracold ^Li and ^'^K atoms e.g. can be tuned into the unitary limit by 
using a magnetic-field Feshbach resonance. Recently measured values for scatter considerably 
and have large error bars: 

§ =0.5 1(4)11, 0.46;i2|3^,0.32+l«@. (8-3) 
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Figure 6: Results for Eq/E^^"^ versus Fermi momentum kf. For comparison we show the results for 
FP 1981 [Hi, APR 1998 [0, CMPR v6 and v8' 2003 [Pi. SP 2005 pi and GC 2007 



Earlier experiments tend to yield larger value for t, indicating the need of further experimental 
studies. 



There have been numerous analytic calculations of (§, see [38] for a recent review. The ob- 
tained values for vary roughly from 0.2 to 0.6. To get a nuclear lattice EFT picture of the neutron 
matter in the unitary regime we simulate the ground state of 8, 12 and 16 neutrons in a box of 
length L = 10, 12 and 14 fm using Monte Carlo. In Fig. |6| we show ground energy ratio Eq/Eq^^ 
in dependence of Fermi momentum kp. For comparison we also show earlier phenomenological 
calculations. Our predictions seem to be consistent with the earlier results. We find a good fit to 
the lattice data using 

Eo/El'^ ~ § - + 0. 16 kp reff - (0.5 Ifm^)^ . (8.4) 

l^F ^scatt 

The results from the fit are 

<^~0.31 and <§i~0.81. (8.5) 



9. N^LO three-body forces 

At N^LO three-body forces start to show up which depend on two constants. We fit these 
LECs from neutron-deuteron scattering data in the spin-1 /2 doublet channel and the triton binding 
energy. To describe the neutron-deuteron scattering we use standard Liischer formula ^2, 33, p4| ]. 



Finite volume spectrum was generated with Lanczos diagonalization method. In Fig. ^ we show the 
S-wave phase-shift in spin- 1/2 doublet channel versus the square of the relative momentum and the 
triton binding energy versus the length of the lattice box. One observes a very natural convergence 
pattern in our simulations with increasing chiral order. Probing the triton binding energy we see 
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Figure 7: In the left panel: neutron-deuteron S-wave scattering phase-shifts in the spin- 1/2 doublet channel 
versus the square of relative momentum. For completeness we show experimental data |44| for proton- 
deuteron and neutron-deuteron scattering. In the right panel: triton binding energy versus the length of the 
lattice box. 



from Fig. ^ that at the box length ~ 15 fm the volume dependence already becomes very small 
and the binding energy approaches its physical value. This is consistent with our expectation that 
the volume dependence in nuclear lattice EFT simulations should become weak for L 20 fm. In 
Fig. ^ we show the S-wave phase-shifts in the spin-3/2 quartet channel versus the square of relative 
momentum. This channel was not taken into account in the fit procedure. Again we observe a very 
nice convergence with increasing chiral order. Our predictions are located between the proton- 
deuteron and neutron-deuteron experimental data. Since the isospin-breaking was not taken into 
account in our simulations the results are very satisfactory. 

As a first Monte-Carlo simulation of N^LO lattice EFT we studied the binding energy of ^He. 
The length of the box was chosen L = 16 fm. In Fig. ^ we show the resulting binding energy of the 
^He system 

{^4\eM-tH/2)Hcxp{-tH/2)\W,) 

^^^"^^ = (^4|exp(-.//)|^4) ^'-'^ 

versus Euclidean time t. Our Monte-Carlo simulations overpredict the physical binding energy 
with subtracted Coulomb-effects by 5%. This is consistent with the expected theoretical accuracy 
of our simulations. 



10. Summary and outlook 

The results of our studies demonstrate that lattice EFT is a promising tool for a quantitative 
description of light nuclei. At leading order binding energies and radii of nuclei up to "^He are 
reproduced with the accuracy 5 . . .30%. At NLO, 9 LECs were fitted to the NN scattering phase- 
shift using the spherical wall method which is best suited to measure phase shifts and mixing 
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Figure 9: Binding energy expectation value of ^He versus Euclidean time t. Plot produced by Monte-Carlo 
simulation with the box length L— 16 fm. 



angles for nonrelativistic point particles on the lattice. With the NLO EFT action, we studied dilute 
neutron matter close to the unitary limit. We performed Monte-Carlo simulation with = 8, 12 
and 16 neutrons in a box of length L = 10, 12 and 14 fm. Our simulation results are consistent with 
earlier phenomenological determinations. We also presented the first analysis of N^LO lattice EFT. 
At this order, the two LECs entering the three-body force were fitted to neutron-deuteron scattering 
data and the triton binding energy. In our first N^LO Monte-Carlo simulation we calculated the 
binding energy of ^He. Our simulations overpredict the physical binding energy of "^He by ~5% 
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which is within the expected accuracy of our lattice simulations. 

In the future, we plan to perform N^LO Monte-Carlo simulations of light nuclei and probe 
neutron matter with larger number of neutrons in a box. 
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